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Recently, the extremely sensitive torsion-rotation transitions in methanol have been used to set a 
tight constraint on a possible variation of the proton-to-electron mass ratio over cosmological time 
scales. In order to improve this constraint, laboratory data of increased accuracy will be required. 
Here, we explore the possibility for performing high-resolution spectroscopy on methanol in a Stark- 
deflected molecular beam. We have calculated the Stark shift of the lower rotational levels in the 
ground torsion-vibrational state of CH3OH and CD3OH molecules, and have used this to simulate 
trajectories through a typical molecular beam resonance setup. Furthermore, we have determined 
the efficiency of non-resonant multi-photon ionization of methanol molecules using a femtosecond 
laser pulse. The described setup is in principle suited to measure microwave transitions in CH[30H[ 
at an accuracy below 10~^. 
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I. INTRODUCTION 



Theories that extend the Standard Model of particle 
physics have presented scenarios that allow for, or even 
predict, spatial-temporal variations of the constants of 
nature 1 . Possible variations of the fine structure con- 
stant, a, representing the strength of the electromagnetic 
force, or the proton-to-electron mass ratio, fi, a mea- 
sure of the strength of the strong force, lead to shifts 
in the spectra of atoms and molecules. Many studies 
have been devoted to observe these shifts. By compar- 
ing metal absorptions in the spectra from distant quasars 
with the corresponding transitions measured in the labo- 
ratory, Webb et al. [2] found evidence that suggests that 
the fine structure constant, a, has a smaller value at high 
redshift. In later work, this variation was interpreted as 
a spatial variation of a 3 . In parallel, laboratory ex- 
periments on earth are used to probe possible variations 
in the current epoch. Compared to their astrophysical 
counterpart, their advantage is their great accuracy, re- 
producibility and unequivocal interpretation. By com- 
paring transitions in different isotopes of dysprosium, a 
possible variation of the fine structure constant was found 
to be < 2.6 X 10~^^/yr !4j. Whereas atomic spectra are 
mostly sensitive to variations in a, molecular spectra can 
be used to detect a possible variation of fi. The most 
stringent independent test of the time variation of /i in 
the current epoch is set by comparing vibrational transi- 
tions in SFg with a cesium fountain, which has resulted 
in a limit for the variation of Afi/fi of 5.6xl0^^^/yr 5j. 
Tests of /z- variation on cosmological time scales have been 
performed by comparing spectra of molecular hydrogen 



measured in the laboratory with those observed at red- 
shifts z — 2 — 3, corresponding to a look-back time of 
10- 12Gyr, constraining A^i/fi < 10~^ [B]. The most 
stringent limit on a variation of /i in the early universe are 
set by Bagdonaite et al. [7] from comparing absorptions 
by methanol in objects at a redshift of 0.89, correspond- 
ing to a look-back time of 7Gyr, with laboratory data. 
The uncertainty in the constraint derived by Bagdonaite 
et al. is dominated by the precision of the astrophysical 
data. However, when more accurate astrophysical data 
become available, the error in the laboratory data will 
become significant. In this paper, we investigate the pos- 
sibilities to increase the precision of selected microwave 
transitions in methanol. We focus on the four transitions 
in CH3OII observed by Bagdonaite et al., and two tran- 
sitions in CD3OII that - provided that the precision is 
significantly enhanced - might be used for a laboratory 
test of the time variation of /i. 

Line centers of methanol transitions in the microwave 
region are typically obtained from absorption mea- 
surements in a gas cell, resulting in (Doppler-limited) 
measurement uncertainties around 50 kHz correspond- 
ing to a relative uncertainty of ~10~'^ [S]. For a lim- 
ited number of lines higher resolution data was ob- 
tained by a pulsed molecular beam Fabry-Perot Fourier- 
transform microwave spectrometer of the Balle-Flygare 
type [3], reaching accuracies around 20 kHz [T(T. Us- 
ing a beam-maser setup, two single methanol transitions 
were recorded with relative accuracies of ^10~^ [TT]. So 
far, this is the only study that was able to (partly) re- 
solve hyper-fine structure in methanol. All these studies 
are based on detecting absorption or emission of the mi- 
crowave field. A significantly higher precision seems fea- 
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sible in a Rabi-type setup using lasers to state-selectively 
detect the methanol molecules. Unfortunately, so far no 
suitable state-selective detection scheme for methanol has 
been demonstrated. The only study that reports the de- 
tection of methanol by resonance-enhanced multi-photon 
ionization (REMPI), involved either the repulsive 3s Ry- 
dberg state or one of the 3p Rydberg state; both resulting 
in broad unresolved bands ^12j . 

Here, we explore the possibility for detecting methanol 
molecules using a femtosecond laser, while relying on 
inhomogeneous electric fields to separate the different 
quantum states present in the beam. This paper is or- 
ganized as follows: In Sec. |ll]we discuss the energy level 
structure of methanol, and review the origin of the large 
sensitivity coefficients that are found in this molecule. 
Furthermore, we outline the procedure that was adopted 
to calculate the Stark interaction for methanol. In 
Sec. |III| we simulate trajectories of methanol molecules 
through a typical beam resonance setup, using the de- 
rived Stark shifts as input. In Sec. IV we present 
measurements that determine the efficiency of ionizing 
methanol molecules using femtosecond laser pulses. Fi- 
nally, in Sec.jVj we use the simulations and measured ion 
yield to estimate the expected accuracy of the described 
beam resonance setup. 



II. THEORY 

A. Torsion-rotation levels in methanol 

Methanol is the simplest representative of the class of 
alcohol molecules and consists of a hydroxyl (OH) group 
attached to a methyl group (CHg). The CO bond that 
connects the two parts of the molecule is flexible, allowing 
the methyl group to rotate with respect to the hydroxyl 
group. This rotation is hindered by a threefold potential 
barrier with minima and maxima that correspond to a 
staggered and eclipsed configuration of the two groups, 
respectively. For the lowest energy levels, the internal ro- 
tation or torsion is classically forbidden and only occurs 
due to quantum mechanical tunneling of the hydrogen 
atoms. In order to account for this additional degree 
of freedom, the normal asymmetric top Hamiltonian has 
to be augmented with a term that describe the torsion 
motion. To simplify the calculation, the coupling be- 
tween overall and internal rotation is partly eliminated 
by applying an axis transformation to the coordinates 
of the Hamiltonian (the so-called " Rho- Axis Method" or 
RAM). In the rho-axis method, the full torsion- rotation 
Hamiltonian for methanol becomes [T51: 



RAM 



(1) 



where iiftors, ^^rot, ^^cd and i?int represent the torsion, 
overall rotation, centrifugal distortion, and higher-order 
torsion-rotation interaction terms, respectively. This 



Hamiltonian is implemented in the BELGI code [13] that 
we have used to calculate the level energies of methanol. 
The current version of the code was modified and im- 
proved by Xu et al. [5] in a number of ways useful for 
treating the large datasets available for the methanol 
molecule. Furthermore, the code has been optimized to 
make it faster and a substantial number of higher order 
parameters has been added. Using a set of 119 molecu- 
lar constants for CHg OH from Ref. [S] and 54 constants 
for CD3OH from Ref. ^14^, the lower energy levels are 
found with an accuracy < 100 kHz. The Hamiltonian of 
Eq. ^ is diagonalized in a two-step process to obtain 
the torsion-rotation energy levels [T5]. In the first step, 
the torsional Hamiltonian (i?tors) is diagonalized: 



Fip^-pJ,y + V{a), 



(2) 



where F is the internal rotation constant, p is the ratio 
of the moment of inertia of the methyl top relative to 
the moments of inertia of the molecule as a whole, and 
V{a) is the internal rotation potential barrier, is the 
internal rotation angular momentum, is the projection 
of the global rotation on the z molecular axis and a is 
the torsional angle. The eigenvalues obtained after this 
first step are the torsional energies for each torsional level 
that are characterized by the quantum numbers K, ut 
and (7 = (A species) or ±1 (E species). The A and 
E symmetry species can be considered as two different 
molecular species in the same sense as ortho- and para 
ammonia. The torsional eigenfunctions can be written as 
linear combinations of the basis wave functions 1151: 
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where k is an integer. In the second step, the rotational 
factor \K) in the eigenfunctions from Eq. Q is replaced 
by the full symmetric top wave function | J KM) to gener- 
ate the basis set used to diagonalize the remaining terms 
of the Hamiltonian from Eq. ([T]), i.e. i?ioti ^^cdi and i/int- 
The lowest energy levels of CHg OH and CD3OH are 
calculated at zero electric field using the molecular con- 
stants from Refs. [H [M], and are shown on the left and 
right-hand side of Fig. [l] respectively. The arrangement 
of energy levels within a symmetry state resembles that 
of a prolate symmetric top, with the difference being 
that every K ladder obtains an additional energy offset 
due to the K dependent tunneling splitting. As a con- 
sequence, certain states in neighboring K ladders may 
become nearly degenerate. It was shown by Jansen et 
al. [T^ and Levshakov et al. [17] that transitions be- 
tween these nearly degenerate states are very sensitive 
to possible variations of the proton-to-electron mass ra- 
tio, /X. The sensitivity of a transition with a frequency v 
is defined as 
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FIG. 1. Level scheme of the lowest torsion-rotation states in the torsion-vibrational ground state (vt = 0) of CH3OH (left 
panel) and CD3OH (right panel). The energies are calculated at zero electric-field strength using the molecular constants from 
Refs. [SI [14] and are given with respect to the zero point of the torsional well. The levels are labeled by Jk (indicated on the 
left-hand side of each level). For the A levels the so-called parity quantum number (-1-/-) is also indicated. Arrows in the panel 
on the left-hand side mark transitions that were used in Bagdonaite et al. [7] to constrain Afj,/fi. The shaded area in the panel 
on the right-hand side highlights the near degeneracies between the 2o E, li E, and 22 E levels in CD3OH. 



The overall rotational energy of the molecule scales with 
its rotational constants and is thus inversely proportional 
to the reduced mass of the system. Therefore, pure rota- 
tional transitions have a sensitivity coefficient, if^ = — 1. 
The torsional energy arises from the tunneling effect and 
- similar to the inversion splitting in ammonia - depends 
exponentially on the effective mass that tunnels. For 
the normal isotopologue of methanol, the sensitivity co- 
efficient of a purely torsional transition has a value of 
= —2.5. Note that, due to symmetry, such transi- 
tions are not allowed in methanol. An interesting effect 
occurs for transitions between different K ladders. In 
this case, part of the overall rotational energy is con- 
verted into internal rotation energy or vice versa. When 
the energies involved are rather similar - i.e., when the 
levels are nearly degenerate - this results in enhanced 
sensitivity coefficients. These enhancements occur gen- 
erally in every internal rotor molecule, but because of a 
number of favorable properties, the effect is exceptionally 
large in methanol [IB] . 

The arrows in the panel on the left-hand side of Fig. [l] 
mark transitions that have been observed in the study of 
Bagdonaite et al. j^. The two transitions in the K — 
ladder are pure rotational transitions and have a sensitiv- 
ity coefficient of i^^ = — 1. The transitions between the 
K — and \K\ ~ 1 ladder have sensitivities of if^ — —7.4 
and —33. The shaded area in the panel on the right-hand 
side of Fig.[l]highlights the near degeneracies present be- 



tween the 2oE, liE, and 22 E levels in CD3OH. The 
22 O li E transition and the li O 2q E have sensitiv- 
ity coefficients of 330 and -42, respectively. Thus, if fi 
increases, the frequency of the 22 O li E' transition be- 
comes larger while the frequency of the li -H> 2o i? tran- 
sition becomes smaller. By comparing these two transi- 
tions over a number of years, a possible variation of fi 
can be constrained or measured. 



B. Stark effect in methanol 

In order to calculate the Stark shift on the energy levels 
of methanol for different values of the electric field, we 
have included the Stark Hamiltonian -ffstark hi^ the second 
diagonalization step of the belgi code. The Stark term 
is given by 

-ffstark = — Me-E, (5) 

where fie is the body-fixed electric dipole moment vector 
and E the electric field. In our calculation, fie is rep- 
resented by the body-fixed dipole moments fj,a and /if, 
along the a and b axes of the RAM frame, respectively. 
For CH3OH, fia = 0.889 D and fib = -1.44 D 0, while 
for CD3OH, fia = 0.8340D and fit = -1.439D [19]. The 
matrix elements of i?stark are taken from Eqs. (l)-(4) of 
Kleiner et al. [IHI ■ Since the Stark effect induces nonzero 
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matrix elements (J|i?stark|>/±1), the Hamiltonian matrix 
was extended to include those interactions. For a given 
J value, the Hamiltonian matrix only has off-diagonal 
blocks involving the nearby J ± 1 states. Each J, K en- 
ergy level is split into 2 J + 1 components, characterized 
by Mj, the projection of the total angular momentum J 
along the laboratory axis Z. Mj is the only good quan- 
tum number in the presence of the electric field, hence 
the basis used for calculating the Stark eS'ect for a certain 
Mj level includes all states with this specific Mj. 

Our approach was tested by comparing the Stark shifts 
calculated by the modified version of belgi with the re- 
sults of a perturbation-like calculation using a code writ- 
ten by Lees and coauthors [21]. For J levels up to J = 2, 
the ratios of the Stark shifts obtained by these two meth- 
ods are close to unity at lOkV/cm and range from 0.7 to 
1 at 50kV/cm. 

In the left and right-hand side of Fig. [2] energies of the 
near-degenerate 3_i E and 2q E levels of CH3OH and 
2qE, liE, and 22 E levels of CD3OH, respectively, are 
plotted as function of the applied electric-field strength. 
The lower panels of Fig. [2] show the corresponding effec- 
tive dipole moments /^cff, defined as 



tark 



d\E\ 



(6) 



where VFgtark is the Stark shift of the quantum state in 
an electric field of magnitude \E\. The effective dipole 
moment is a measure for the orientation of the molecule 
in a specific state, i.e., the expectation value of the dipole 
moment in a space-fixed axes system. Note that in the 
Hamiltonian used, \J,+K,—Mj),E levels are degener- 
ate with I J, -K, +Mj),E levels (if K and Mj ^ 0) 
and, as these levels have different parity, this degener- 
acy results in a effective dipole moment that is non-zero 
at OkV/cm. This non-physical result disappears when 
high-order couplings are incorporated [Hj. In Table |l] 
the effective dipole moments of a selection of levels are 
listed. 



III. TRAJECTORY SIMULATIONS 

In this section, we use the obtained effective dipole mo- 
ments to simulate molecular trajectories through a typi- 
cal beam resonance setup [S^. We consider a molecular 
beam apparatus that consists of (i) a collimation section; 
(ii) an electric defiection field for state preparation; (iii) 
a microwave cavity; (iv) a second electric deflector for 
state-selection and (v) a detector. In the calculations, 
we assume a methanol beam with a mean velocity of 
800 m/s. This beam is collimated by two 0.6 mm diame- 
ter skimmers that are separated by 500 mm. The deflec- 
tion fields used for pre and post state selection consist of 
a cylindrical electrode and a parabolically shaped elec- 
trode to which a voltage difference is applied. The elec- 
tric field inside such a defiector was analyzed by de Nijs 



and Bethlem in terms of a multipole expansion It 
was shown that the most optimal defiection field is cre- 
ated by choosing a field that only contains a dipole and 
quadrupole term, the strength of which is represented 
by ai and 03, respectively. In our calculations, we take 
ai=18 kV and 03=8.6 kV while all other expansion terms 
are set to 0. Furthermore, tq that characterizes the size 
of the electrodes is taken as 3mm (i.e., the distance be- 
tween the electrodes is 2ro = 6 mm) . This results in an 
electric field magnitude of 60 kV/cm on the molecular 
beam axis. For CH3OH the deflectors are assumed to 
be 200 mm long, while for CD3OH a length of 150 mm is 
used. Both deflectors are oriented in the same direction, 
i.e., a molecule that is deflected upwards in the first elec- 
tric field, will again be deflected upwards in the second 
electric fleld, if it has not undergone a microwave transi- 
tion. The force on a molecule inside one of the deflection 
fields is found from: 



-McffVl^l, 



(7) 



where the gradient of the electric field is calculated from 
an analytical expression derived in de Nijs and Beth- 
lem [21]. The two defiection fields are separated by 
200 mm. After the last deflection field, the molecules 
travel 200 mm further before being ionized in the focus 
of a femtosecond laser. The total flight path adds up 
to 1350 mm for CH3OH and 1250 mm for CD3OH. The 
200 mm free flight between the deflection fields can be 
used to drive a microwave transition. 

In our simulations, the trajectories of typically 5 x 10^ 
molecules are calculated. The position and velocity 
spread are sampled randomly from a Gaussian distribu- 
tion. The initial quantum state is sampled from a Boltz- 
mann distribution that assumes a temperature of either 
5 or 10 K. As mentioned earlier, A and E type methanol 
should be considered as two different molecular species. 
In A-type methanol, the total spin is / = |, whereas 
in E'-type methanol I — hence, the nuclear spin de- 
generacy for A-type methanol is twice that of i?-type 
methanol. However, this is matched by the presence of 
both K > and K < levels of E symmetry [25]. In 
our calculations, we assume equal numbers of molecules 
of A and E species for CH3OH, whereas the A/E ratio 
in CD3OH is taken to be 11/16 

Fig. [3] shows the distribution of the effective dipole 
moment of the molecules. The black curve (a) shows 
the distribution when the deflection fields are off and the 
laser focus is situated at the molecular beam axis. The 
red (b) curve shows this distribution when the deflec- 
tion fields are on and the laser focus is situated 2.2 mm 
above the molecular beam axis. It is observed that only 
molecules are detected in states that have an effective 
dipole moment of — 0.5 x 10~^ cm~^/(kV/cm). With 
these settings, the resolution of the selector, A/icg, is 
0.2x 10"^ cm~^/(kV/cm). In these simulations, the Mj 
state of the molecule is assumed to be preserved - i.e.. 
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FIG. 2. Calculated level energies and effective dipole moments of the 3_i E and 2o E levels in CH3OH (left panel) and 2o E, 
li E, and 22 -E levels in CD3OH (right panel) as function of the applied electric field. 



it is assumed that a homogeneous electric field is ap- 
plied between the deflectors to keep the molecules ori- 
ented. However, in our experiment the region where the 
microwave excitation takes place should be completely 
shielded from external magnetic or electric fields. In this 
case the different Mj levels are degenerate and, in the 
worst case scenario, the Mj distribution after the inter- 
action zone is randomized completely. The results of a 
simulation that assumes complete de-orientation is shown 
as the blue curve (c) in Fig. [sj Clearly, Mj scrambling 
greatly reduces the effectiveness of the state selection. 
Note that due to hyperfine splittings, the degeneracy of 
the different levels at zero electric field is lifted. The hy- 
perfine structure has not been taken into account in our 
simulations as it is largely unknown, in fact, resolving 
the hyperfine splittings is an important motivation for 
this study. 

Fig. |4] shows a simulation of what would be observed if 
the vertical (y) position of the laser focus is scanned, 
i.e., the spatial distribution of the methanol beam in 
the direction of field-gradient. The red curve shows the 
distribution that is observed when the microwave field 
drives 50% of the molecules from the Jk = 3-i E to 
Jk = 2o i? level, and vice versa. We assume that the mi- 
crowave field only drives AMj = transitions. The black 
curve shows the distribution when no micro- wave field is 



present. In both cases the Mj distribution after the in- 
teraction zone is assumed to be randomized completely. 
Both curves are normalized to the number of molecules 
that are detected at y = when the defiection fields are 
turned off. We define the count rate with the microwave 
field being on or off resonance as Ron and i?off, respec- 
tively. The count rate of the undefiected beam is defined 
as Rq. In order to observe if a transition has occurred, 
we look for a difference between i?off and i?on, i-e., the 
difference between the black and red curves in Fig. |4] 
This difference is shown in the lower panel of Fig. [4j If 
we choose the vertical position of the laser focus to be 
2.3 mm above the molecular beam axis, Ros/Ro= 0.062 
and i?on/-Ro=0.048. This results in a normalized aver- 
aged count rate, R = ^{RoS + i?on)/-Ro — 0.06, and a 
difference in count rate, AR = {Ros — Ron)/Ro = 0.013. 
In Sec. JVj it is shown that the accuracy scales with 
\AR\/VR. For the 3_i — >■ 2o transition, this number is 
equal to 0.056. Table|n]hsts AR and \AR\/^/R for the se- 
lected transitions assuming a rotational temperature of 5 
and 10 K. The errors in AR and \ AR\/V~R are estimated 
to be 10-20% for CH3OH and 20-40% for CD3OH. For 
the selected transitions in CH3OH \ AR\/VR is 0.03-0.08 
at 5 K and slightly less at 10 K. For the two considered 
transitions in CD3OH, these numbers are even less favor- 
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TABLE I. Effective dipole moments of selected torsion-rotation states of CH3OH and CD3OH. 
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able, which is shghtly surprising given that the effective 
dipole moments of the levels involved are rather different; 
the different Mj states of the Jk = 22 -E level all have a 
positive Stark shift, the Mj states of the Jk = "^o E level 
have a negative Stark effect, whereas the Mj states of 
the Jk ~ li E level have virtually no Stark shift. How- 
ever, the fact that the effective dipoles of the different 
Mj states within the levels varies considerably, combined 
with Mj scrambling in the interrogation zone, compli- 
cates state-selection. Furthermore, the population of the 
levels involved is rather small. 



IV. NON-RESONANT IONIZATION OF 
METHANOL USING FEMTOSECOND LASER 
PULSES 

In order to estimate the efficiency of non-resonant 
multi-photon ionization of methanol by a femtosecond 
laser, we have performed test measurements in an exist- 
ing molecular beam machine described elsewhere |27j . A 



supersonic molecular beam is prepared by expanding a 
mixture of ~0.5% CH3OH in argon through a 200 /im di- 
ameter pulsed piezo nozzle [55] operating at a repetition 
frequency of 10-1000 Hz. The pressure behind the nozzle 
was kept below 1 Bar in order to prevent cluster forma- 
tion. The molecular beam passes a 1.5 mm and a 1.0 mm 
skimmer, before entering the detection region where it is 
intersected at right angle with a focused (/ = 500 mm) 
femtosecond laser beam. The total flight path from the 
nozzle to the detection zone adds up to 160 mm. The 
produced ions are accelerated towards a position sensi- 
tive microchannel-plate detector mounted in front of a 
fast phosphor screen where they are counted. The fem- 
tosecond pulses are generated by a commercial amplified 
regen laser system (Spectra Physics Spitfire) that pro- 
duces a 1 kHz pulse train and is tunable around 800 nm 
with an output energy of about 800 /xJ and a duration 
of 120 fs. The fundamental 800 nm light is doubled in 
a BBO crystal to obtain a 400 nm laser pulse with an 
energy of about 100 /iJ. As the ionization potential of 
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TABLE II. Simulated normalized difference in count rate, AT? = {RoB — Ron)/Ro, and \AR\/V~R with R = |(-Roff + Ron)/Ro, 
for selected transitions in methanol at a temperature of 5 and 10 K. The fourth and seventh columns indicate the y position 
of the focused laser beam used that results in the maximum ratio of AR/y/R. The third column lists the sensitivity of the 
transition for a variation of the proton-to-electron-mass ratio, K^. 
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FIG. 4. Upper panel: simulated spatial distribution of the 
methanol beam along the y axis when the microwave field 
is off resonance (black curve) and on resonance with the 
3-1 2o E transition (red curve), assuming a rotational tem- 
perature of 5 K. Both curves are normalized to the number of 
molecules that are detected at y = when the deflection fields 
are turned off. Lower panel: difference between the count rate 
on and off resonance. 



FIG. 3. Distribution of the effective dipole moments of the 
methanol molecules in a molecular beam with a temperature 
of 5K. The black curve (a) shows the distribution when the 
deflection flelds are off and the laser focus is situated at the 
molecular beam axis. The red (b) and blue (c) curve show 
this distribution when the deflection flelds are on and the laser 
focus is situated 2.2 mm above the molecular beam axis. For 
the red curve, Mj is assumed to be preserved in the fleld-free 
region, whereas for the blue curve complete Mj scrambling is 
assumed. The curves are given an offset for clarity. 



methanol is 10.84 eV p5^, at least four photons of S.lOeV 
are required to ionize it. Table [III] hsts the number of de- 
tected parent ions per second using an intensity of 20 
or 100 /xJ per pulse at a IkHz repetition rate. Although 
the ion yield can probably be increased by further opti- 
mahzation of the molecular beam parameters and focal 
properties of the laser beam, it seems unlikely that this 
changes the obtained results by more than a factor of 
2. The fact that the nonresonant multi-photon ioniza- 
tion rate does not follow a /" dependence, with / the 
laser intensity and n the number of photons absorbed. 
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can be understood by considering the intensity profile 
of the laser focus. At a certain laser intensity, /gat, the 
ionization rate at the beam waist reaches unity, and the 
signal increases due to the increase in volume for which 
I > -^sat- For comparison, the ion yield for fluoromethane 



CHgF is also listed in Table 



III Note that the ionization 
and at least 



potential of fluoromethane is 12.50 eV 
five photons are required to ionize it. 



TABLE III. Number of methanol and fluoromethane parent 
ions resulting from a 400 nm 120 fs laser pulse at a repeti- 
tion rate of 1 kHz. For this measurement, methanol and fluo- 
romethane molecules were seeded in argon with relative con- 
centrations of 0.5 and 5%, respectively. 



molecule 




R (ions/sec) 






20 J /pulse 100 M.I /pulse 


CH3OH 


10.84 


1200 20 000 


CH3F 


12.50 


3500 25 000 



^ Data from Lias et al. [29] 



As expected, the Allan deviation becomes infinite when 
the count rates on and off resonance are equal, while it 
scales with 1 / \/Ros when i?on is zero (background free) . 
The extra factor of arises because we spend half of 
the time on signal and half of the time off signal. The 
values for \/\R\/\fR are given in Table i?o, the count 
rate when the deflection fields are turned off, can be es- 
timated from the test measurements. Assuming that the 
density in the molecular beam decreases quadratically 
with the distance from the nozzle, the expected signal at 
a distance of 1350 mm behind the nozzle is about a factor 
of 64 smaller than that obtained in the test setup, hence 
i?o = 3.1 X 10^ ions/sec. For the 3_i E to2QE transition 
at 12 GHz, the Allan deviation is 7.0x 10"''/Vr. In order 
to achieve a fractional accuracy of 10"^ a measurement 
time of about 80 minutes would be required. Similar re- 
sults are found for the other transition in CII3OII listed 
in Table [llj The expected accuracy for CD3OII is less 
favorable. 



V. ESTIMATED ACCURACY 

From the measured ion yield and the results from the 
simulations, we can now estimate the accuracy that can 
be obtained in the described beam machine. In general, 
the accuracy of a frequency measurement depends on the 
Q factor of the system and the signal to noise ratio, S/N . 
In an experiment that relies on counting individual ions 
the accuracy is expressed by the Allan deviation [31)] 

where defines the duration of one measurement cycle, 
T is the total measurement time, v is the frequency of 
the measured transition and lS.v the width of the spectral 
line. In our simulations, we assume that the microwave 
cavity (or microwave cavities, if we use a Ramsey type 
setup) has a length of 160 mm while the molecular beam 
has a velocity of 800 m/s. This implies that the total 
interrogation time for the measurement is 200 /is, corre- 
sponding to a spectral width of ^4 kHz. For the 3_i E 
to 2o £^ transition at 12 GHz, this results in a Q-factor 
of 3 X 10^. The S/N depends on the number of ions 
that are detected per cycle. If we spend half the dura- 
tion of a measurement cycle on and off resonance, the 
number of ions that contribute to the signal is given by 
S — Tc\Rofi — i?on|/2. Assuming Poissonian statistics, 
the noise on the total number of detected ions is given 
by = ^rc(i?ofr + i?on)/2. Thus the Allan deviation 
becomes; 

""^ Q |i?off- i?o„| V r " Q |Ai?| 



VI. CONCLUSION 

In this paper we estimate the feasibility of performing 
high resolution microwave spectroscopy on selected tran- 
sitions in methanol using a Rabi-type molecular beam 
setup in combination with a femtosecond laser. We have 
adapted the belgi programme to calculate the Stark ef- 
fect for the different isotopologues of methanol. The cal- 
culated Stark shifts are reasonably large, typically be- 
tween -1 and Icm^^ in a field of lOOkV/cm for the dif- 
ferent rotational states. Thus, the molecules can be easily 
manipulated using modest sized deflection fields. Due to 
the small rotational constants of methanol, many states 
are populated even at the low temperatures that can be 
obtained in supersonic beams. With the resolution ob- 
tained in a typical molecular beam defiection setup, it is 
not possible to select individual quantum states. Typi- 
cally, 5 or more states are present within the laser focus. 
Furthermore, in the field free region used for inducing 
the microwave transition, the Mj quantum number is 
not preserved, leading to a reduction of the state purity. 
Note that, Mj-scrambling is both a blessing and a cure, 
on one hand it reduces the observed difference between 
Ron and i?off ; on the other hand it ensures that all hyper- 
fine levels are addressed without he need to change the 
position of the laser focus. Mj-scrambling can be en- 
forced by rapidly switching the second defiection field on 
entrance of the molecules '31]. From simulations, we find 
that the differences in signal on or off resonance, AR, 
are in the range of 0.01 — 0.02 for the selected transi- 
tions in CH3OH, and slightly less for the two considered 
transitions in CD3OH. In order to estimate the detec- 
tion efficiency of methanol molecules using a femtosec- 
ond laser, we have performed test measurements in an 
existing molecular beam machine. Using a laser power 
of 100 /iJ per pulse at 1 kHz repetition rate, 2 x 10* ions 
per second were detected at a distance of 160 mm behind 
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the nozzle. From these numbers, the described molecular 
beam deflection setup seems suited to measure microwave 
transitions in CH3OH at an accuracy below 10"^, al- 
though the required measurement times are rather long. 
Currently the laser system is being upgraded to provide 
peak intensities that are 5 times higher than used in this 
study. This should result in an increase of the ion yield of 
at least a factor of 25, making this experiment perfectly 
feasible. Due to the small differences in count rate, how- 
ever, it might be tedious to find the optimal position of 
the laser focus. For a laboratory test using CD3OH an 
accuracy well below 10^^*^ would be required. This seems 
unlikely to be achieved in a beam machine as considered 



here. 
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